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Abstract 

A non perturbative numerical method for determining the discrete spectra is de- 
duced from the classical analogue of the Schrodinger's equation. The energy eigenvalues 
coincide with the bifurcation parameters for the classical orbits. 

1 Introduction 

One of known curiosities of the one-dimensional Schrodinger's equation is its reduciblity 
to the 1-st order differential equation of Riccati. Given the eigenvalue problem for the 
Schrodinger's wave packet ip(x): 

-\^4>(x) + [V(x)-E]4>(x) = (1) 

(where m = fi = 1), the formal substitution: 

^(x) = e^ x) (2) 
leads to the classical Riccati equation for f(x) = <f>'(x): 

f(x) + f( x y = 2[V(x)-E] (3) 

Despite its apparent simplicity, (3) is one of non-trivial and persistently returning prob- 
lems in mathematical physics. In fact, even the field equations of the General Relativity 
might be viewed as a 'spiritual analogue' of the Riccati eq. (see e.g. [1]). An obvious idea 
would be to replace completely the spectral problem (1) by its 1-st order equivalent (3). 
The difficulties, though, are formidable. In the first place, the ansatz (2) is impractical if 
ip(x) is a real sign changing function. Moreover, (3) has a tendency to create singularities 
of f(x), which appear even for V(x) = const, and correspond to the nodal points of ip(x) 



1 



= Vo =>■ f{ x ) = ~k tan(kx), where k = \/2(E — Vq)]. The subtler ansatz of 
Aharonov and Au [2] liquidates only a part of this difficulty. 

A different chapter in these attempts was open by the angular analog of the Riccati eq.(3) 
found by Drukarev [3] and Franchetti [4]. Yet, their works appeared in wrong time (wrong 
place?) and did not receive the attention which they deserved. Drukarev and Franchetti 
themselves, apparently, have formulated their equations just to facilitate the calculation of 
the phase shifts (which they did, but the achievement later faded together with the whole 
"phase shift trend"). As a result, the multiple forms of the Riccati eq. remain a kind 
of sophisticated curiousm, basically known but scarsely applied. Notable exceptions are 
the techniques of evaluating the number of energy levels for one dimensional Schrodinger's 
operator developed in [4, 5, 6, 7, 8], and specially works of Calogero [9, 10]. Yet, the 
capacities of the method to determine the energy spectra are far from explored. Our purpose 
here is to complete the story, by studying the exact numerical consequencies of the angular 
equations [3, 4]. We shall show, that their numerical integration is, most likely, one of the 
most powerful approximate methods to determine the eigenvalues of one-dimensional (or 
spherical) spectral problems, which might reduce the traditional perturbation calculus, at 
least for 1-dimensional eigenvalue problems, into a kind of museal relict. 



2 The classical model of the Schrodinger's equation. 



Curiously, the meaning of our algorithm is best seen by forgetting completely about the 
quantum mechanical sense of the Schrodinger's equation and sticking to its classical equiva- 
lent. This point of view, though seldom applied, has some notable traditions (see, e.g. [11]). 
Its most provocative expression was the description of the Saturn rings as the band spectrum 
of the Schrodinger's operator [12]. 

To perceive the classical sense of (1) denote the variable x by t and call it time; put also 
q= ip(t), p= ip'(t). Equation (1) then becomes: 



2[V(t) - E]q 



(4) 



o?q dp 
dt ^' dt 

Note, that (4) coincide with the canonical equations for the pair of classical variables q,p 
defined by a time dependent Hamiltonian: 



H(t) 



P 



+ [E- V(t)W 



(5) 



representing the classical oscillator with a time dependent quadratic potential V(q, t) 
[E — y(t)]q 2 . The classical phase space trajectory of (4-5): 



q(*) 



q(t) 
P(*) 



(t G R) 



(6) 



'paints' a detailed image of the Schrodinger's wave function ip(x) and its first derivative ip'(x). 
This fact has been decisive to represent the spectral bands of the Schrodinger's operator as 
the stability bands of a classical system (the Saturn's rings! [12]). It will be as decisive in 
our description of point spectra. 
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3 The bifurcations. 



Assume for simplicity, that none of E, V(t) is positive, and V(t) vanishes outside of a 
finite interval [a, b] £ R (V(t) = for t < a and t > b). The equations (4-6) thus describe a 
classical point moving under the influence of a constant repulsive potential -Eq 2 , corrected by 
an attractive term — V(t)q 2 . For E < 0, the typical trajectory tends to ±oo as t — > ±oo (for 
great \t\ the repulsive force dominates, expulsing the classical point to infinity). Occasionally, 
however, a curious dynamical phenomenon occurs: some trajectories, emerging from the 
phase space origin (q=p=0) at t = — oo, by a rare coincidence, acquire a momentum sufficient 
exactly to return asimptotically to the origin, against the repulsive forces. This phenomenon, 
extremally unstable, represents precisely the eigenvectors of (f) [ip(x) — > as x — > ±oo], 
i.e. the most stable forms of quantum motion. The effect, in a sense, looks as a classical 
skill game (whose "goal" is to collocate the material point at the repulsion center; see [13]). 
To explain the mechanism of the game, some geometry elements on the phase plane V are 
needed. 

Due to the linearity of (4-5) the phase vector q(t) at any t £ R depends linearly on the 
initial vector q(a) for any a £ R, i.e.: 

q(t) = u(t,a)q(a), (7) 

where the 2x2 simplectic real matrix u(t, a) is called the transference or evolution matrix. 
The eqs. (4) traduce themselves easily into the first order matrix equation for u(t,a): 



du 
~dt 



A(t)u(t, a) 



(8) 



where: 



A(t) 



(9) 



1 

2[V{t)-E] 

Mathematically, the matrix equation (8-9), is neither easier nor more difficult to solve 
than the original Schrodinger's equation (1). Yet, it leads to geometric pictures which 
facilitate the solution of the spectral problem (1). 

Consider first of all the region outside of the potential well, = R\(a, 6), where V(t) = 0. 
The generator A(t) in is constant: 



A(t) = A 

leading to the explicit solutions: 

q(*) 



Note, that A fulfills: 



1 
2 |E| 



; A (*- a )q(a) 

; A(f-&) q ( 6 ) 



for t < a or t > 6, 



for t < a 
for t > b 



(10) 



111 



A 



2 IE! 1 



12) 



and so, it has two eigenvalues \± = ±\/2 \ E\ and eigenvectors e±: 
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eigenvalues 



eigenvectors 



A+ = +a/2TE 



+ a/2 \E\ 



(13) 



A_ = -y/2\E\ 



2 E 



Henceforth, for t G (i.e. in absence of V(t)) the motion on the phase plane V preserves 
the directions e± producing a continuous squeezing: the direction e + expands while e_ 
exponentially shrinks as f -> +oo (inversely for t — > — oo). (The squeezing is typically 
generated by the hamiltonians of repulsive oscillators, while the attractive ones generate the 
phase space rotations). 

The phase trajectory (6-7), in general, diverges in both +oo and — oo. However, excep- 
tions exist. If the initial phase vector q(a) is proportional to e + , then in agreement with 
(11) q(t) vanishes as exp[(t — a) y/2 \E\] for t — > — oo. In turn, if q(&) is proportional to e_, 
then q(t) vanishes as exp[—(t — b) y/2 \E\] for t — > +oo. A number E < is an eigenvalue of 
the Schrodinger's operator in (1) iff there exists a non-trivial trajectory vanishing on both 
extremes t — > — oo and t — > +oo. This can happen if and only if the evolution between t = a 
and t = b brings e + into a vector proportional to e_, i.e: 

u(b, a)e + = Const x e_ (14) 

To monitor the phenomenon, consider the integral trajectories of (4-5) for different values 
of E < [or alternatively, for different amplitudes of Let 7 + , 7_, — 7 + , — 7_ be the 

radial |-lines sticking from the phase space origin in the directions of e + ,e_,— e + , — e_ and 
let V + and V~ denote two open connected components of V\I- containing 7 + and — 7 + 
respectively (see Fig.l). In absence of V(t), V + and V~ are evolution invariant; as t — > +oo, 
all the trajectories in V + escape to infinity tending asimptotically to 7 + -axis while all the 
trajectories in V~ tend to — 7 + . Consider now the trajectory which departs from q( — oo) =0 
crossing q(a) = e + (E). If V = 0, q(t) = exp[(t — a) \J2 |75|]e + (i.e., our trajectory escapes to 
infinity exactly along the 7 + -axis). If |V(t)| is non zero but small in [a, &], the escape along 
7_|_ (due to the squeezing/expansion) is corrected by — V(t)q 2 which generates a clock-wise 
rotation in V . Yet, once the rotation dies out at t > 6, the squeezing takes over and drives 
the trajectory back again to 7 + as t — > +oo. If |V(t)| is stronger [or alternatively, if \E\ is 
weaker] the trajectory will deviate more from 7 + and return to it slower (Fig.l). Finally, 
for sufficiently small \E\ (or sufficiently strong the rotation caused by V(t) will bring 

q(&) to 7_ (see Fig.l). Instead of returning to 7 + , the phase point will fall down right to the 
origin along the (shrinking) direction 7_, drawing a (normalized) eigenvector of (1). For \E\ 
still smaller, q(&) will cross to V~ and the trajectory will alter its asimptotic behaviour: it 
will tend now to — 7 + as t — > +oo. For \E\ further decreasing [or |V| increasing], q(&) will 
circulate on V crossing subsequently — 7_, 7_, — 7_, 7_, . . .. Every such event will produce 
a closed trajectory (an eigenvector!), coinciding with an abrupt change of the asymptotic 
type. The eigenvalues of the energy operator are thus interpretable as bifurcation values (i.e. 
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the values of E for which the classical orbits produce bifurcations) [An idea arises, that the 
discrete spectrum of (1) could be simply defined in terms of bifurcations. Such definition 
would not require linearity, and so, could be easily extended to the non-linear Schrodinger's 
operators (see the discussions in [14, 15] ]. 




Figure 1. The metamorphosis of the classical trajectory (4-5) for varying E 
and fixed V(t) (qualitative picture). As E raises to zero from below, the defor- 
mation due to the rotating term — V(t)q 2 expands clockwisely around the phase 
space origin, crossing several times the "shrinking axis" e_(£"). At each new 
intersection a bifurcation occurs, producing a new closed orbit interpretable as 
an eigenvector of the Schrodinger's equation (1). The trajectory transforma- 
tions are pictured in the moving frame of the 'squeezing axis' and represent as 
well the bifurcations which must occur for a fixed E < and variable V(t). 

What not less interesting, the bifurcation values can be determined by a rather simple 
algorithm, independent of the traditional perturbation methods. 

4 The angular Riccati equations. 

Since the vector norms in (14) are irrelevant, (14) can be conveniently written in terms of 
an angular coordinate. Indeed, introduce: 
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q = / ocosa, p = / osina (15) 

The canonical equations become: 

p cos a — ap sin a = psina (16) 
p sin a + ap cos a = 2[V(t) — cos a (17) 

where p and a mean the time derivatives. Curiously, the equation for the angular variable 
separates. Multiplying (16) by — sin a, (17) by cos a (or viceversa) and adding one gets the 
1-st order differential equation for a alone: 

a = 2[V(t)-E]cos 2 a-sin 2 a (18) 

and simultaneously: 

p/p= [V(t)~ E + l/2]sm2a (19) 

The equation (18) was first derived by Drukarev [3] and Franchetti [4] (by using purely 
algebraic arguments). Its link to the original Ricatti idea is immediate. In fact, (18) implies: 

d 

— t&na + t&n 2 a = 2[V(t)-E] (20) 

The advantages of the Drukarev- Franchetti form over the traditional Riccati equation 
(3), however, are that: 1) the eq. (18) can be solved for arbitrary a's without leading to 
singularities; 2) it offers a clear geometric idea of the spectral condition, and 3) it leads to 
an easy numerical algorithm. Indeed, notice, that the phase vectors e± define the angles: 



®±(E) = ± arctan \J2 \E\ (21) 



6(E) = a + (E) - a_(E) = 2 arctan ^2 \E\ (22) 

In agreement with (11), the trajectory tends to zero in — oo iff a(a) = a + + kir } and it 
tends to zero in +oo iff aib) = a_ + rmr. The trajectory tends to zero for both t — > ±oo 
iff the time evolution (18) converts the initial angle a(a) = a + into aib) = a_ ± n7r. 
Thus, to check whether a number E < is an energy eigenvalue it is sufficient to solve the 
differential equation (18) with E fixed and with the initial condition a(a) = a + (E) } finding 
aib) = aib, E). Whenever the total angular change Aa = a(a) — aib) in [a, b] fulfills: 

Aa = 6(E) + nvr, n = 0,1,2,... (23) 

E belongs to the discrete spectrum. A convenient form of (23) involves the spectral defect 
angle defined as a difference between the desired angle a -(E) and the achieved angle a(b } E): 

T(E) = a_(E) - a(b } E) = Aa- 6(E) (24) 
The condition (23) then tells: 

Y(E)=mr n = 0,1,2,... (25) 
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An immediate generalization of the conditions (23-25) is obtained for V(t) constant but 
not necessarily vanishing outside of (a, b): 

™={$; tin ^ 

The motion (4-6) has then two different constant generators A(a) and A(&) for t < a and 
t > b and the formulae (23-25) hold after substituting — E\ or \V(b) — E\ instead of 

\E\ in the expression (21) for a + and a_ respectively. 

Observe, that (15) is not the only way to separate the angular part of (4). The "P-plane 
is a simplectic space without a natural measure of distances and angles. Henceforth, instead 
of introducing the polar variable (15) straighforwardly, one might as well introduce new 
canonical coordinates q',p', and only after define the angular variable on q',p' plane. Some 
profits of this freedom were explored by Calogero [9]; here let us notice only the plausible 
form of the spectral condition if the angular variable is introduced by: 

q = A/) cos a , , 

p = A -1 /} sin a 

where A = (2 | ^E 1 1 ) ^ . The motion equations then read: 

a = ^2 \E\ cos 2a + ^2/ \E\ V(t) cos 2 a (28) 

p/p = {^2\E\ + V{t)y/\E\ /2 ) sin 2a (29) 
and the limiting angles are independent of E (a± = ±^). The spectral condition becomes: 

Aa = (n + i)7r (30) 

in a visible reconciliation between the oscillator spectrum and the Sommerfeld quantization 
conditions. 



5 The numerical algorithm. 

The importance of the angular variables for the structure of the bound states was detected 
already in papers devoted to the phase shifts [17]. The defect angle equivalent to (24) was 
introduced indeed by Calogero, and used to evaluate the number of the energy levels for the 
1-dimensional and radial Schrodinger's eigenproblem [10]. Yet, the efficiency of the method 
to determine the exact eigenvalues somehow escaped attention (perhaps, due to a general 
fascination by the perturbative methods!). The main reason, why (18) is so easily applicable, 
is: 

Theorem 1. For any bounded, piece-wise continuous V(t) with a compact support in 
[a, &], T(E) is a strictly increasing function of E for E < V(a). 

The proof is deduced from two observations. (1) If ait) and a'it) are two solutions of 
(18) with a(a) < a' (a) then aib) < a'ib) [Indeed, otherwise there would be a point c £ [a, b] 
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with o(c) < o'(c) contradicting the uniqueness of the solutions of (18)]; (2) If a\{t) and a 2 (t) 
are solutions of two equations of form (18) with two different parameters E = E\ and E = E 2 
respectively, then Oi(a) = 0:2(0), an d E\ > E 2 (and) 01(a) < o 2 (a) [the proof involves only 
the elementary analysis of the Cauchy eq.(18)]. The observations (1) and (2) imply now that 
a(b } E) for an angular trajectory starting in o(a) = a + (E) is a monotonically decreasing 
function of E, and the proof is completed by noticing that o_(_E) is increasing. 




Figure 2. The defect angle T(E) = a_ — a(a,E) for two cases of truncated 
oscillators: (a) b = —a = 2 and (b) b = —a = 4. The intersections of the 
"stepping" functions T(E) with the lines n x 1 give the eigenvalues of the 
Schrodinger problem. 



As an illustration, we have used (25) to determine the energy levels for the truncated 
1-dimensional oscillator potential: 



\uj 2 x 2 for \x\ < a 



V(x)={ 2 (31) 



^u 2 a 2 for \x\ > a 
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The limiting angles are a± = ± arctan \J w 2 a 2 — 2E. We have determined the angular 
function T(E), < E < V(a), for w = 1, a = 2, and for a = 4, by integrating numerically 
the angular equation (18) (see Fig. 2). It yields the 2 energy levels for the oscillator truncated 
at a = 2, and 8 energy levels for the oscillator truncated at a = 4, all calculated with accuracy 
up to 10~ 10 . Curiously, the obtained eigenvalues are very close to the first 2 and 8 levels 
of the exact oscillator, respectively, E n = n + | (indeed, even the last eigenvalues of the 
truncated potential (31) differ rather little from the orthodox E\ = 1.5, and E7 = 7.5). 
Comparing to Ritz method, a basic advantage of our algorithm is its essential simplicity (no 
need to waste skills inventing an adequate class of the test functions). A notable advantage 
is that the spectral function T(E) is unstable and changes very quickly when crossing the 
sequence of critical values T = mr (n = 1,2, . . .) (see Fig. 2). Thus, even a very little error 
in E in the vecinity of an eigenvalue, traduces itself into a visible effect in T, significantly 
improving the accuracy (compare Calogero [16, p 274]). This instability is caused by the 
fact that the energy eigenvalues correspond to the bifurcations of the orbits and the final 
point aib) deflects very fast when E crosses the bifurcation value. In limit b = +00, a(b, E) 
would be discontinuous and T(E) would be an exact step function! 



6 The "! eigenproblem" . 

Our method, till now, concerns only the potentials constant outside of finite intervals (lim- 
ited, non-singular wells). Could it tell something about more general V{t)1 Consider any 
continuous V : R — > R such that: 



V± = liminfV(t) > -00 (32) 

In the traditional approach to the spectral problem, the main effort is to find the non- 
trivial solutions of (1) vanishing on both extremes t — > ±00 (which exists only as an excep- 
tion!) Following the observations of Sec. 5 we propose to reduce the solution to two minor 
steps, each one interpretable as "| spectral problem": 1) for any E find the special solutions 
of (1) which vanish for t — > — 00 (the "left eigenvectors"); 2) find the solutions of (1) which 
vanish for t — > +00 (the "right eigenvectors"). 

While the difficulty of solving the complete spectral problem is formidable, one seldom 
pays attention to the fact, that every "half of it" has a solution always, defining some left 
(right) decaying branches for (1) for any E < V- (E < V + ). For big \t\ they provide the 
asimptotic cues e±(t, E) and asimptotic angles a±(t, E) adequate to replace the fixed vectors 
e±(E) (13) and angles ct±(E) in the algorithm of Sec. 5. Indeed, one has: 

Lemma Let f3 :R— >R be a continuous real function with: 



liminf /3(t) > r/ 2 > (33) 

t^ + OO 

Then the 2-dimensional solution space S of the 2nd-order differential equation: 

" S = mq(t) (34) 
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has a 1-dimensional subspace ^_ of solutions which vanish for t — > +00 and are square 
integrable in [0,+oo). 

Proof. Let N be a number such that fiit) > rj 2 for t > N. For t > N, the material point 
q(t) moves under the influence of the repulsive ellastic force 

q > r/ 2 q (35) 

Consider now an integral trajectory of (33) which satisfies the initial condition: q(iV) = 
C > 0, q(N) = j]C > 0. Using (35) one easily shows that q(t) is positive, monotonically 
increasing and: 

q(t) > q(N)e^- N ^ = Ke nt (36) 
The method of variation of constant then provides a new, linearly independent solution: 

, , f +co dr ^ f +co dr f +co dr e""* , , 

q - (t) ^ (t) i ^ 2q(i) i M=L W)-ct (37) 

which spans the desired subspace 
An immediate consequence is: 

Theorem 2. Let V(t) be a continuous potential in the 1-dimensional Schrodinger's 
equation (1) and assume (32) holds. Then for every E < V_, the 2-dimensional solution space 
of (1) contains a 1-dimensional subspace S + (£') of solutions vanishing for t — > — 00, square 
integrable in ( — 00, 0], while for any E < V + it contains a 1-dimensional subspace S_(£') of 
solutions vanishing for t — > +00, square integrable in [0,+oo). A number E <min(V_, V+) 
is a point of the discrete spectrum of (1) if the subspaces S_(£') and S + (£') coincide. 

Proof is obtained by applying the Lemma to t and —t with fiit) = V(t) — E and 
< J] 2 < V± — E respectively. 

In the Appendix we have collected the asimptotic forms of both "decaying branches" of (1) 
for several V(t) including the oscillator and Coulomb potentials (see Appendix). These cues 
(the "|-eigenvectors"), permit to solve, with any desired accuracy, the traditional spectral 
problem for a class of infinite potential wells. Let V(t) be one of these potentials and let 
Eq < Voo. The energy eigenvalues E < E (if they exist), can be then determined by the 
following algorithm, generalizing the one of Sec. 5: 

i) one fixes a finite interval [a, b] CR such that V(t) — E > n > and the asimptotic 
cue expressions be valid for t < a and t < b; 

ii) one uses the "vanishing cues" q±(t,E) to define two special angles: 

ct-(a,E) = arctan[p_(a, E)/q_(a } E)] (38) 
a+(b,E) = nrctcin[p + (b,E)/q + (b,E)] (39) 
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characteristic for the trajectories vanishing at t — > — oo and t — > +00. 

iii) one integrates the equation (18) in [a, b] with the initial condition a(a) = 
a + (a } E) finding aib) = a(b } E) and determining the defect angle: 

T(E) = a_(b } E)-a(b } E) (40) 

iv) whenever: 

T(E)=mr, re = 0,1,2,... (41) 
the number E < E belongs to the point spectrum of the Schrodinger's operator. 

The non-trivial part of the method, of course, is to find the vanishing cues [18]. However, 
as the interval [a, b] can be arbitrarily wide, it is enough to know the asimptotic expressions. 
Now, if (32) holds and E < E 0} the cues are monotonic, without zeros for t < a and t > b; 
one can thus use the anzatz (2-3). Note that while the general solution / of the Riccati 
equation (3) depends on one arbitrary constant, and (tipically) diverges as < -> ±00, the 
fit) of vanishing cue has no such arbitrariness, and can be determined with any desired 
accuracy by applying the known iterative or asymptotic methods (see Appendix). For the 
oscillator and Coulomb potentials the cues ("^-eigenvectors) are already known in form of 
the asimptotic series for the confluent hypergeometric equation, vanishing either at t — > + or 
t — > 00 (however, we prefer to represent them in the form q(t) = exp{f(t)}, since fit) = tan a 
defines the asimptotic angles). 

What quite essential, once the cues q± (and the angles a±) are determined, they can 
be used not only to find the spectrum for one particular potential V(t) but simultaneously, 
for an entire class of potentials which share the asimptotical behaviour of V(t) (and can be 
arbitrarily deformed in any finite region). Moreover, given the "left cues" for one potential 
V\(t) and the "right cues" for another potential V^t) the method can be as well used to 
determine the spectrum of any V(t) sharing the asymptotic behaviour of V\{t) for t — > — 00 
and of V2(t) for t — > +00. 

For curiosity, we have used the asimptotic angles calculated in our Appendix to determine 
the spectrum for the "hybrid oscillator" (Fig. 3) not so easily treatable by either perturbative 
or variational methods: 

{\x 2 for x > 
(42) 
\x 2 for x < 

Note, that our method permits to solve as easily any other hybrid or deformed cases (as 
e.g. an oscillator affected by an arbitrarily high potential barriere in the middle, etc.) 
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E2=l . 667399152170 



E1=0. 996678885237 



E0=0. 339309660629 



Figure 3. The hybrid oscillator potential. The dashed lines represent the 
standard energy eigenvalues for each separated oscillator, and the solid lines 
represent the actual eigenvalues for this potential. 



7 Singular and radial wells. 

The physically important wells not only extend to infinity, but can have singularities in the 
finite region. The typical case is the 1-dimensional eqs. (1) obtained after the separation of 
the angular variables in the Schrodinger's eq. in R 3 with a radial potential <f>(r). By denoting 
?/>(x) = R(r)Y(6 } (f>) and assuming Y(6,<f)) to be an eigenfunction of the square angular 
momentum L 2 , one ends up with the 1-dimensional eigenvalue problem for uit) = tRit): 



d 2 



_ + [2V (t) + 



2E]t 







(43) 



dt 2 1 v y t 2 

whose eigenvectors are the trajectories vanishing for t — > 0+ and square integrable in any 
[6, +oo) (b > 0). Similarly as before, each of the asimptotical conditions, typically, can be 
satisfied by solutions of (1) for any E. The existence of the right cues (solutions vanishing 
for t — > +oo) is assured by Theorem 2. In turn, the existence of the zero-cues (solutions 
vanishing at 0+) is the consequence of the following elementary theorem: 



Theorem 3. Let cf>it) be a continuous real function in (0,+oo), satisfying: 

^ (f)(t) > -k(t) 
in a certain subinterval (0,(5 > 0), where k : (0,(5) — »R fulfills: 

'■5 ft 



k(t) > and 



o 



k(t')dt' < + 



oo 



(44) 



(45) 



Then the 2-dimensional space ^ of the functions u : (0, +oo) — >R which solve the 2nd-order 
differential equation: 
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d?u 

- — + <f>(t)u = (46) 

must contain a 1-dimensional subspace S of solutions which vanish for t — > + . 

Proof. We shall stick to the classical image of Sec. 2 putting uit) =q and uit) =p; (46) 
then paints the motion of a classical point mass m = 1 under the influence of the ellastic 
force — 2cf>{t)(\{t) . Choose now a £ (0, 6) and consider the integral trajectories of (46) defined 
by the initial condition p(a) = q(a) = (q(a) arbitrary). The assumption (45) permits one 
to evaluate the lower bound of q(t) in (0,a). Indeed, 

m > -k(t)q(t)^ (47) 

pa pa pa 

=>• q(a) - q(t) = / q(t')dt' > - / k(t')q(t')dt' > -q(a) / k(t')dt' => (48) 
Jt Jt Jt 

=> q(t) < q(a) J k(t')dt' (49) 

a pa 



where 



Henceforth, if a is small enough, so that K(a) < 1, then q(t) has a positive lower bound in 
(0,a). If now <f)(t) < in (0,a), then q(t) < q(t) > = q(a) q(t) monotonically 
increases in (0, a) tending to its positive lower bound at t — > + . In turn, if cf>it) > in (0, a), 
then q(t) < and q(t) monotonically decreases in (0, a) tending to its (finite or infinite) 
upper bound in (0,a). In both cases, the standard method of the "variation of constant" 
provides the new solution of (46): 

* (i) = q(i) l^ = i '" (i) ui w ) (53) 

where hit) = l/q(t) is positive, monotonic and bounded in (0,a). Note that the expression 
(1/t) J Q h{r)dr is the average value of h(r) 2 in (0,t). If <f>(t) is positive, and h(t) increasing 
(<^q(t) decreasing) the fraction in the square parenthesis in (53) is < 1; but if cf>it) is negative 
and hit) decreasing, it tends to 1 from above as t — > + . In both cases, for t small enough, 
there is a constant N > such that qo(t) < Nt, and so, qo(t) defines the desired subspace 
S of solutions vanishing as t — > + .n 

In the Appendix we have collected examples of asimptotic cues at + for some tipical 
singularities. 

To check the results, we have used the right cues it — > +00) and the + cues of the 
Appendix to determine the discrete spectrum for the Coulomb well V(t) = — j, obtaining 
the first 10 levels of the hydrogen atom with accuracy to 10~ 10 : 



q(a)-q(t)= / q(t')dt' < q(a)/ / k(t n )dt n dt' < K(a)q(a) (50) 

Jt Jt' 

^q(t) > [1-A»]q(a), (51) 

pa pa 

K(a) = k(t')dt'dt. (52) 
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E 


= -0.5000000000 


Ei 


= -0.1250000000 


E 2 


= -0.0555555555 


E 3 


= -0.0312499999 


E 4 


= -0.0200000000 


E 5 


= -0.0138888888 


E 6 


= -0.0102040816 


E 7 


= -0.0078125000 


E 8 


= -0.0061728395 


E 9 


= -0.0050000000 



Once the asimptotic expressions were cross-examined, we have used the already tested 
right cues of the oscillator and the + cues of the Coulomb singularity to find the principal 
series of energy levels for the "hybrid well" : 

V(r) = -- + ~w 2 r 2 (54) 
r 2 

sometimes considered as a candidate to describe the quark confinement [19]. Taking w 2 = 
0.00005 we could observe, that a very weak oscillator potential cancels the condensation 
of the hydrogen energy levels for E — > 0~, providing a continuous transition to an equally 
spaced spectrum for E > —0.01 (see Fig. 4). Let us also notice that a close relative of our 
method has been succesfully used to study the logarithmic wells [20]. 




Figure 4. A hybrid well: the kind which could resemble the quark confinement. 

We are tempted to predict that in not so far future, the perturbative methods of de- 
termining spectra (at least for the 1-dimensional Schrodinger's operator) will be almost 
forgotten and a part of pacience nowadays devoted to special functions (permitting to solve 
only exceptional problems!) will be invested into building up a bank of data about the 
asimptotic behaviours and "vanishing cues". Once these data will be precise enough, the 
task of determining the spectra of arbitrary potentials of known asimptotic types will become 
a question for the pocket calculators. 
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Appendix: The evaluation of the vanishing cues 

Below, we look straightforwardly for / = tan a defined by the Riccati equation (3) and 
yielding q(t) = ef ^ dt which vanish either for t — > ±oo or t — > + . 



A Harmonic oscillator 

The solution to the Riccati equation 

f + f 2 = wH 2 - 2E 
which yields q(t) — > for t — > ±oo, can be found in form of an infinite series: 

/(*) = -wt + a ^~ i 



where 



Explicitely: 



a — a 2s — 0, 



8 = 



E 1 
w ~ 2 : 



a s+1 = ^- j-(> - l)a s _! • "> | 



(s > 2) 



fit) 



-wt + 
1 



~E 


1' 


1 1 

1 + 2w 


~E 


r 




~E 


3" 


w 


~ 2_ 


w 


~ 2_ 




w 


~ 2_ 



1 



2w 2 



E 1 
w ~ 2 



E 3 

w ~ 2 



E 

2 

w 



The asimptotic form of the "^-eigenvectors" for t — > ±oo is: 



-2('#-aU-4f#-f)(^-2] 



The consistency with the textbooks expressions for ip(t) in terms of the confluent hypergeo- 
metric function is easily verified. 
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B The Coulomb potential 



The series solution lor 



[yielding qo(t) — > when t — > 0+] turns out to be 



Explicitely 



8 = 



= 111 _ _J_ _ 2^(/ + l) 2 + l 
/o " t /+1 (2/ + 3)(/ + l) 2 



Henceforth, the "|-eigenvector" vanishing for t — > 0+ is: 

q ~ e i+1 2(2 i+3 )( i+ i)2 I +- 

consistently with the well known expressions in terms of the confluent hypergeometric series. 
For t — > +00, by similar arguments, / has the form 



f /oipi 1 1 1 2l(l + l)\E\ + y/2\E\-l l 



yielding the "|-eigenvector" 

q_ ~ f?W e h(i/t) 

where has the form of an anlytic series vanishing at £ = 0. However, we have found 
that if the integration interval is wide enough, a simpler asimptotic expression gives good 
results: 

/ 2 /(/+ 1) 

f- = -\ -2E--+ { ^ } 



t ' t 2 

and the only problem we have is to integrate the angular equation (18) in a sufficiently long 
range, to have the cue without nodal points. The integration can be simplified by changing 
variables from t to x = t / (1 -\- 1). 



C The Yukawa potential 

This case is very similar to the preceding one, though no longer treatable in terms of hyper- 
geometric series. The Riccati equation for fit) = tana(t) reads 
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Our Theorem 3 assures the existence of the solution of (1) vanishing at zero, of the form 
ipo(t) = ex p{fo(t)} where: 

i=0 

After calculations 

_ _-]_ _ 2(A - E)(l + l) 2 - 1 

00 ~ 7 + 1' 01 ~ (2/ + 3)(/ + l) 2 

■1 I 2(-A) s+1 



2/ + S + 3 I (s + 
Henceforth, the vanishing cue qo(t) is: 



U>1) 



r=0 



2(X-E)(l + lf-l 2 _ 

q (t) ~ t l+1 e 1+1 2 ( 2 i+ 3) (i+ i)2 



(no longer representable in terms of the confluent hypergeometric series!} 
For t — > +00, we have used a finite approximant: 



f- 



1(1+1) 2 w 

V ; -2E- -e" At 
t 



changing simultaneously the integration variable to x = t/(l + t). 



D A hybrid "quark potential" 

The Riccati equation (3) for the effective potential 



V(i) = -\+ l -wH 

yields the following fo(t) for t — > + : 

l + l 1 2£(/ + l) 2 + l 



2.2 



fo 



2t 2 



2E(l + l) 2 + l f2 + 



t l + l (2/ + 3)(/ + l) 2 (2/ + 3)(/ + l) 3 (/ + 2) 



1 



21 + 5 
The vanishing cue 



w 



AE(l + l) 2 + 2 



(2£(/ + l) 2 + l) 



(2/ + 3)(/ + l) 4 (/ + 2) (2/ + 3) 2 (/+l) 4 



t 6 + . . . 



q = e -«TT^W 

where has the form of an analytical series. 

Meanwhile, the solution vanishing at t — > +00 corresponds to f-(t) in the form 



/_ = -wt + 
The vanishing cue is 



E 1 
w ~ 2 



111 1 

t wt 2 2w 



q_ ~ e 2 

where is an analytic series vanishing at £ = 



~E 


1 




"£ 


3" 


w 


~ 2_ 




u; 


~ 2_ 


(-- 

V w 


-§) e fc(i/*) 
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